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Введение. Статья посвящена разработке расчетной модели 
бесконечного радиального подшипника, смазываемого рас- 
плавом легкоплавкого покрытия, с учетом зависимости вязко- 
сти смазочного материала от давления. Рассмотрены парамет- 
ры, обусловленные расплавом легкоплавкого покрытия под- 
шипниковой втулки от удельной теплоты плавления, и зави- 
симости вязкости смазочного материала от давления. Дана 
оценка влияния указанных факторов на основные рабочие 
характеристики трибосистемы. Цель работы — формирование 
уточненных расчетных моделей бесконечных радиальных 
подшипников, работающих в режиме гидродинамического 
смазывания при наличии смазочного материала и расплава 
легкоплавкого покрытия подшипниковой втулки, с учетом 
зависимости смазочного материала от давления для примене- 
ния в инженерной практике. 

Материалы и методы. Предложены новые математические 
модели, описывающие движение несжимаемого смазочного 
материала в приближении для тонкого слоя, уравнение нераз- 
рывности и выражения скорости диссипации энергии для 
определения профиля расплавленной поверхности легкоплав- 
кого покрытия подшипниковой втулки с учетом влияния ряда 
дополнительных факторов. Выполнен сравнительный анализ 
вновь полученных и уже имеющихся результатов, что под- 
твердило приближенность новой модели к реальной практике. 
Результаты исследования. Разработаны новые многопара- 
метрические выражения для основных рабочих характеристик 
рассматриваемой пары трения, учитывающих зависимость 
смазочного материала от давления при наличии смазочного 
материала и расплава легкоплавкого покрытия подшипнико- 
вой втулки. Дана оценка влияния параметров, учитывающих 
целую гамму переменных факторов, обусловленных распла- 
вом поверхности легкоплавкого покрытия подшипниковой 
втулки от удельной теплоты плавления. 

Обсуждение и заключения. В предлагаемой работе обобщено 
влияние пока еще не исследованных факторов, что суще- 
ственно усложняет задачу, но делает ее решение универсаль- 
ным и востребованным в современных трибоузлах. Результа- 
ты численного анализа полученных теоретических исследова- 
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ний показывают, что опоры скольжения, работающие на рас- 
плаве легкоплавкого покрытия, обладают аномально низким 
коэффициентом трения (зависимость коэффициента трения от 
параметра, обусловленного расплавом, близкая линейной). 
Полученные результаты могут быть использованы в условиях, 
когда подача смазочного материала связана с трудностями, в 
частности, в таких отраслях, как машиностроение, авиастрое- 
ние, приборостроение ит. д. 


Ключевые слова: гидродинамика, радиальный подшипник, 
вязкий несжимаемый жидкий смазочный материал, расплав- 
ленная поверхность подшипниковой втулки, зависимость 
вязкости смазочного материала от давления. 
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Введение. В современном машиностроении трибоузлы для новых машин проектируются, как правило, с уче- 
том увеличения статических и ударных нагрузок, действующих на опоры скольжения. Следует отметить, что одним из 
важнейших конструктивных элементов подшипников жидкостного трения является смазочная среда. 

Одним из методов решения конструктивно-эксплуатационных задач может быть применение смазывания рас- 
плавом легкоплавкого покрытия подшипниковых втулок. 

В настоящее время существует значительное число технологий диффузионной металлизации, большинство из 
которых являются сложными, требуют применения дорогостоящего оборудования, что делает их промышленно нере- 
ализуемыми. На этом фоне выгодно отличается технология диффузионной металлизации сталей из среды легкоплав- 
ких жидкометаллических растворов. Ее перспективность связана с простотой реализации, возможностью использова- 
ния стандартного термического оборудования, совмещения данной технологии с термической обработкой покрывае- 
мых изделий и т. п. Однако самое главное — возможность получения качественных изделий и регулирования свойств 
покрытий за счет варьирования технологических режимов. 

Диффузионная металлизация из среды легкоплавких жидкометаллических растворов основана на явлении се- 
лективного изотермического переноса металлических элементов на металлические материалы [1—4]. Технологически 
данный процесс осуществляется путем погружения и выдержки изделий в легкоплавком расплаве, содержащем в рас- 
творенном состоянии элемент, на базе которого формируются диффузионные покрытия. 

Использование данной технологии позволяет получать на поверхности изделий однокомпонентные и много- 
компонентные диффузионные покрытия на базе различных металлических элементов. 

Смазывание жидкими металлами применяют при температурах, при которых обычные смазочные среды пре- 
терпевают необратимые физико-химические изменения. Преимущество смазывания расплавом состоит в том, что сма- 
зочный материал образуется в области контакта, где это необходимо. Плавление доставляет достаточное количество 
смазочного материала в зону трения, механические и конструктивные сложности, связанные с его подачей, отсут- 
ствуют. Смазывание расплавом изучалось во многих прикладных задачах, в частности в процессах формоизменения и 
резания металлов [5—7]. Большое количество работ [8—14] посвящено гидродинамическому расчету радиальных под- 
шипников бесконечной длины в условиях отсутствия смазочного вещества и учета зависимости вязкости смазочного 
материала от давления. Существенным недостатком рассматриваемой пары трения, работающей на смазывании рас- 
плавом, является низкая несущая способность. Кроме того, процесс смазывания не является самоподдерживающимся. 

Таким образом, разработка расчетной модели подшипников скольжения, работающих на смазочных материа- 
лах в виде металлических расплавов, с учетом вышеуказанных аспектов функционирования, представляет собой одно 
из перспективных направлений теоретических исследований современной трибологии. Последнее определяет новизну 
и актуальность полученного решения. 


Научная новизна предлагаемого решения и уточнения расчетной модели заключается в оценке влияния па- 
раметра, обусловленного расплавом легкоплавкого покрытия подшипниковой втулки, и параметра, характеризующего 
зависимость вязкости смазочного материала от давления, обеспечивающих опорам скольжения аномально низкий ко- 
эффициент трения. 
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Цель работы: формирование уточненной расчетной модели опоры скольжения в режиме гидродинамическо- 
го смазывания при наличии смазочного материала и расплава легкоплавкого покрытия подшипниковой втулки с уче- 
том зависимости вязкости смазочного материала от давления. 

Постановка задачи. Рассматривается модель установившегося движения вязкого несжимаемого смазочного 
материала в рабочем зазоре бесконечного радиального подшипника скольжения, покрытого легкоплавким расплавом. 

Вал вращается с угловой скоростью &), а подшипниковая втулка неподвижна. Предполагается, что простран- 
ство между эксцентрично расположенными валом и подшипником полностью заполнено смазочным материалом, а 
подшипниковая втулка выполнена из материала с низкой температурой плавления. 

Рассматриваются условия, когда все тепло, выделяющееся в смазочной пленке, идет на плавление поверхно- 
сти материала подшипниковой втулки. 

Задается следующая зависимость вязкости смазочного материала от давления: 


Г _ ар’ 
и ше , (1) 
где ио — характерная вязкость; и’ — коэффициент динамической вязкости смазочного материала; р’ — гидродина- 


мическое Давление в смазочном слое; 9 — постоянная. 
Исходные уравнения и граничные условия. Движение смазочного материала описывается уравнением те- 
чения вязкой несжимаемой ЖИДКОСТИ В приближении для тонкого слоя и уравнением неразрывности: 





й ‚ 
у а ду, У, 10 
г и РР 90, (2) 
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Здесь у,,у„ — компоненты вектора скорости смазочной среды; р’ — гидродинамическое давление в смазоч- 
ном слое; и’ — коэффициент динамической вязкости. 


Рассмотрим полярную систему координат с полюсом в центре подшипниковой втулки (рис. 1). 











Рис. 1. Расчетная схема радиального подшипника с легкоплавким покрытием 
Ето. [. Апа?ас тоае! о гай Беаптиз уий ие соайпе 


В этом случае уравнения контура вала расплавленной поверхности легкоплавкого покрытия подшипниковой 
втулки и поверхности подшипниковой втулки, покрытой металлическим расплавом, запишем в виде: 


м =и,(1+Н), "=к, м =к+^.1 (0), (3) 
Г... е :. 
где Н =5с050—— Е? 1? 0+... ==—; И — радиус вала; И — радиус подшипниковой втулки, покрытой металли- 
о 
ческим расплавом; е — эксцентриситет; & — относительный эксцентриситет; }/{ (0) — ограниченная функция при 


9О> [о = 2% | подлежит определению. 


Граничные условия в рассматриваемом случае с точностью до членов О(=2) запишем в виде: 


м =0, у,„=0 при "=н+А (0); 





у =КО, у,=-Оезт0, при г’=к+есо$0; (4) 
! ! Р 
р'(0)=р’(2п)= > 
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Для определения функции /,'/ (0) ‚ обусловленной расплавленной поверхностью подшипника, воспользуемся 


формулой скорости диссипаций энергий: 





Аал’Г(09)х +есоз0 (р 2 
(6) ого [ м (5) 
49 п++2.Г(6) {07 
где Г —— удельная теплота плавления на единицу объема. 
Переход к безразмерным переменным реализуется на основе следующих формул: 
| ‚ шов 
И=н-6, З=Л-Ж; %=Оду у =Оби рр’ р’ =НО: 
’ (6) 


| а 
И =; @=-. 
р 


Выполняя подстановку (6) в систему дифференциальных уравнений (2) и (5), а также в граничные условия (4), 
придем к следующей системе дифференциальных уравнений: 


р. 
ди де 40 ОГ 00 
ал/"(9 1-1 с0$0 - 
ы а — ( 
40 —п, (6) Г Го 


и граничным условиям: 
у=1, и=-1$00, при г=1-1со$0; 








р. (9) 
у=0, и=0 при г= 0-Ф(6); р(0) = Р(2п) = 
РВ 
е ^' 
где =>; =>; Ф(9)=п./ (6). 
Пусть 2 =е. Продифференцировав обе части равенства, получим: 
4 _ ор и, сор _ 142. 
40 40 40 © 40 
Тогда уравнения (7) и (8) примут следующий вид: 
2 
ИВ А (10) 
ди" © 40 О’ 00 
4Ф(9 1-1 5050 (0 = 
р Га (11) 
49 0 [4 
с соответствующими граничными условиями: 
у=| и=-1$10, при г=|-тсо$; 
ар (12) 


у=0, и=0 при г=0-Ф(0); 2(0)=2(2*)=е Р. 
В качестве малого параметра принимая К, обусловленный расплавом и скоростью диссипации энергии, будем 


искать функцию Ф (09) ввиде: 
Ф(0) =-КФ (0)-К”Ф,(0)-К?Ф.(0)-...=Н (13) 
гдеН =-КФ (0)-К* Ф,(0)-К?Ф.(0)-.... 
Граничные условия для безразмерных компонентов скорости и и у на контуре у = -Ф(0) можно записать в 


виде: 


[А д”у 2 
и(0-Н(9)) =у(0) -— —- Н(0)- = Н7(9)-...=0 
Г /ь=0 0т г=0 


2 
и(0— Н(@)) = и(0)- [>= Н(9)- Е Н?(0)-..=0 (14) 
г=0 и=0 


Асимптотическое решение системы дифференциальных уравнений (10)—(11) с учетом граничных условий (12) 
и (14) будем искать в виде рядов по степеням малого параметра К: 
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У=У (7,0) + Ку (г,0) + К’, (г,0)+...; 
и=и (г, 9) + Ки (т,0) + Ки, (х,0)+...; 
Ф( 0) =-КФ, (0)-К°Ф, (9) -К°Ф, (0)-...; 
Ф(0) =-КФ (0)-К*Ф,(0)-К?Ф.(0)-...; 
2=27,+КА (0)+К"7, (0)+К?7. (6)... (15) 
Выполняя подстановку (15) в систему дифференциальных уравнений (10)—(11) с учетом граничных условий 


(12), получим следующие уравнения. 
Для нулевого приближения: 


2 
7 
Об _ 14% 0% 0% _/ 








= - 16 
ди" © 40 00 0 (16) 
с граничными условиями 
и=Ь щ=-П5Ш0 при г=1-1со$6; 
У 0. Ио =0 при г = 0; (17) 
Э 
25 (0)= 2 (2п)=е ?Р. 
Для первого приближения: 
У _ 14. дм _,, 
д? 040’ 00 бб" “ 18) 
аФ. (9 1-1с050 [0 2 
2%] 
40 0 [бд 
с граничными условиями 
ОУ ди 
у =| | .„Ф (0): и =| —® .Ф (0): 
| |. ке * [5 ы 
у=0; м=0 при г=1|-1со$6; (19) 


2 (0)=2,(2^)=0; КФ, (0)=Ка, Ф(0)=Ф(2”)=а. 
Точное автомодельное решение. Точное автомодельное решение задачи для нулевого приближения будем 
искать в виде: 


ОГ 


_ Оуо 


У +Й, (*,0); щш = - +0 (х,0); 


Г 


Чо (7,8) = 9025 81 ооб (20) 


И, (",0)=$(=); (И, (и,0) = (&)-1'(0); 
47 С С 
ба + — ‚ 1(0)=1-псо$6. 
49 й (0) й (0) 
Подставляя (20) в систему дифференциальных уравнений (16) с учетом граничных условий (17), получим сле- 
дующую систему дифференциальных уравнений: 


ЕС» Иась (<) 8) 50 (21) 








и граничные условия: 


$, (0)=0, 9'()=0, &()=—пзте, #(1)=0: 


1 
д (0)=0, $(0)=ь [%(8)4Е=0. (22) 
0 
Непосредственным интегрированием получим: 
— С _ . С > _ 
и (Е) =—2(8' 8), все [а С 6. (23) 
Ра 


Из условия 25 (0) = 2 (2п)=е Р получим следующее выражение: 
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т 49 | | 

® 0 12 (0) 12 1-12 

Сб ао о | Г 
й° (0) 


Определение гидродинамического давления. С учетом (24) для о получим: 


= = : бат 9 1 ти 
и бес О а й = Па о. (25) 
° 172 (8) +654 ) а Я 


Общеизвестно, что в случае решения задачи плоской гидродинамической теории смазки давление определя- 
ется с точностью до произвольной постоянной. Значение этой постоянной С, устанавливается из условия отсутствия 


отрицательных значений давления в смазочном слое. 
Для нахождения значения 2 для первого приближения необходимо сначала определить функцию Ф, (0) | 


Для определения Ф, (0) с учетом уравнения (22) придем к следующему уравнению: 








2 

аФ (6) _ #(0) [45 (8) $ (2), 
= 26 
40 7. [ое | ИВ 

Интегрируя уравнение (26), получим: 
°’^490 9 /^4,а0 9 ^.а0 

0 1 р 3 

= | УОМУХНОМУЯДО ео 
1 т . _ РВ М 
где о (5). в 45-19 |. 5-7 (8)4Е= 6 ” Аз =1(5(5)) 46 =4 


: (28) 
зир р — зир И, Е 
[0-2] ° [0-2л] (2+1? )(1- пс0$ 0) (1-03 0) 64 
Решая уравнения (27)—(28) с учетом Ф (0) = а, ‚ получим: 
4(1+2'] 6бт(21? +1) зт@ 6[1-п?* ]зт 0 
9, (9) =. ыыы ОО а т та |- не. За ЕО. (29) 
(2-1 и (2+1?) (1-1с0$6) (2+1?) (1-1с0$6) 
Точное автомодельное решение задачи для первого приближения будем искать в виде: 
д д 
У ЕРЕ+И (7,0); и! = = +01 (*, 0); 


и 


ци! (7,0) =, (Е); Ре. 
й (1,0) = (5); И (,0) =-й (&)-1" (6); 


- 30 
р а т в. ‚ 1(0)=1-псо$6. Г’ 
40 и (0) 1п?(60) 





Подставляя (30) в систему дифференциальных уравнений (18) с учетом граничных условий (19), получим сле- 
дующую систему дифференциальных уравнений: 


$'(5)=С,; 9"(8), =С; 9 (8)-8Я (&)=0 1 
и граничные условия: 
$, (0)=0, \4(1)=0, #(1)=0, 9(1)=0; 

1 
а (0)=0, 9#(0)=м, [9(2)48=0. (32) 

0 

Непосредственным интегрированием получим: 

а © ня [С 

(=) =—2 (2-5). п(8) 4+ = М Е +М, С =6М. (33) 
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Из условия 7, (0) = 7, (2п) =0 получим следующее выражение: 


_ __ 2 
- _12м (1 п?) 











‚ | (34) 
2+1 
_ ду _ —2 6(1- т?) 
и - р Е в (6) _ в |-1с0$09 (2+1?) @-песз6 
где 
64 (2+1?) 1-м 2 (2+1?) (1-1с0$0) (2+1?) (1—исоз0)° 


С учетом (34) для (1 получим: 
в -«[5 1 (в)-+С, () - лы (35) 


(2-1? )(1-ис0$0 —1с0$0) 


Тогда для Д = о + КИ! получим следующее выражение: 

















Ра 
| 1 5-й: 
И.И — (1+КМ)-+е . (36) 
(2+1? )(1-пс0$0) (1-ис0$0) 
Ра 
а. 
или е “= до (1+ КМ )+е о 
бзш 0 ] 
где Л = ИЕ ИИ (Е ИН 
(2 }(1- 1056) (1-ис0$0) 
и. 
Применяя разложение в ряд Тейлора для функций е ®, е Р , получим: 
2 2 2 
ор 108 |8 | = (1+КМ). (37) 
2 р 2\Р 
р. 3 
Решая уравнение (37) с точностью до членов О(а)}?, О(&), о] 5 | для гидродинамического давления, полу- 
чим 
р р. а?[р.\ 
р=-*-4(1+КМ) о ый В В (38) 
Р р 2\Р 
Результаты исследований и их обсуждение. Перейдем к определению основных рабочих характеристик 
подшипника. 


С учетом (16), (18) и (38) для составляющей вектора поддерживающей силы и силы трения получим: 
3 2 2 
ОЗ 2 Р броблул(1+ КМ (п 1-п +1) 2 
К и о Ст 1 Е ^ Ра. 
0 


бд | 52 (2+2 } (1-12) р р. 














иоОх? 2*| бу, ду 4пиоОг? | -(1+2КМ) 31-1?) (1+ км о? р? 
ий | Е Е рр и | 
0 г=0 г=0 1-п (2-м?) (1-7) 
3 2л Р 
К. т [ Г кывшо-о (39) 
0 0 Р 


Для проверочных расчетов на основе полученных теоретических моделей использованы следующие значения: 
шо = 0,00595 Нс/м?; = 0,3...1 м; = 0,019995...0,0493 м; 
О = 100...1800 с"; 6 = 0,05 : 10 3...0,07 : 10-3; К= 0,00052...0,0000022; 
р Ма: О =55.25 561 Н/м?; М = 0,16...25,6. 


По результатам численных расчетов построены графики, приведенные на рис. 2—5. 
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Рис. 2. Зависимость компонента поддерживающей силы (А, ) от параметра а, 


характеризующего зависимость вязкости от давления, и от ея ‚ характеризующего давление на торцах интервала, при К = 0 
Р 


Ето. 2. ерепаепсе о{зиррогйпе }отсе сотропеип! ( К, ) оп рагатеег а ищсй спагачет2ех ргеззиге-у15со5йу гапо, 


апа оп Е дехстбтз рге5зиге аё йе 1тегуа[ епау а К = 0 


ааа 
АмМАУЩщЮЮ 


``. % 


0,002 0,001 


00040.03 





ме 0,8 01005 


Рис. 3. Зависимость компонентов поддерживающей силы (А, ) от параметра а, характеризующего зависимость вязкости 
от давления, и от К, обусловленного расплавом и скоростью диссипации энергии 
Е!о. 3. Оереп4епсе о? зиррогап® Югсе сотропеп$ (А, ) оп рагатаег а, ус свагасеттез ргеззиге-у1$созИу гайо, 


ап оп К сачзе4 Бу тей ап4 епегоу 4135 1райоп ге 





0,0003 
0,0005 0,000 
к 


Рис. 4. Зависимость коэффициента силы трения от параметра К, обусловленного расплавом и скоростью диссипации энергии, 
и от а, характеризующего зависимость вязкости от давления 


Ето. 4. Оерепаепсе о р7тсНоп югсе асю! оп рагатеег К итсй саизеа Бу тей ап4 епегзу @5я1рапоп гие, 
апа оп а спагаметлтие ргез5иге-у15со5Йу гапо 
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Рис. 5. Зависимость коэффициента силы трения от параметра а, характеризующего зависимость вязкости от давления, 
и от гидродинамического давления р, при К = 0 


Ето. 5. Оерепаепсе о 7споп рюгсе дасюг оп рагатеег о спагаметилив ргеззиге-у5со5Пу гано, апа оп пуагодупатис ргеззиге р аК = 0 


Выводы. Анализ полученных расчетных моделей и графиков позволяет сделать ряд выводов. 

1. Получена уточненная расчетная модель радиального подшипника скольжения, работающего в условиях 
гидродинамического смазывания расплавом легкоплавкого покрытия с учетом зависимости вязкости от давления. 

2. Показан значительный вклад конструктивного параметра К, обусловленного расплавом. С увеличением 
конструктивного параметра К коэффициент трения уменьшается на 60 %, а несущая способность увеличивается на 
16 %. 

3. Зависимость коэффициента трения от конструктивного параметра К, обусловленного расплавом, близка К 
линейной в принятом нами приближении в пределах 0,0014—0,003. 
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